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The reaction pd — * pX at 305 MeV is studied with the aim to search for supernarrow dibaryons, 
the decay of which into two nucleons is forbidden by the Pauli exclusion principle. The experiment 
was carried out at the Moscow Meson Factory using the spectrometer TAMS, which detected the 
scattered proton and another charged particle (either p or d) from the decay of X. Narrow peaks in 
missing mass spectra have been observed at 1905 and 1924 MeV. Comparison of the obtained data 
| with theoretical predictions leads to the conclusion that the peak found at 1905 MeV most likely 

corresponds to a supernarrow dibaryon with the isotopic spin equal to 1. The possible origin of the 
£SJ ' peak at 1924 MeV is also discussed. 

PACS number(s): 14.20.Pt, 12.39. Mk, 13.75.Cs 
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(N ■ I. INTRODUCTION 

The possibility of the existence of multiquark states was predicted by QCD inspired models These works 

C-T) , initiated a lot of experimental searches for six-quark states (dibaryons). Usually one looked for dibaryons in the NN 
channel. In the present work we will consider supernarrow dibaryons (SNDs), the decay of which into two nucleons is 
forbidden by the Pauli exclusion principle ■ In the NN channel such dibaryon states correspond to even singlets 
and odd triplets at the isotopic spin T = 0, and odd singlets and even triplets at T = 1. These dibaryons with the 
mass M < 2mN + Tn v {mjsiijn^) is the nucleon (pion) mass) can decay into two nucleons, mainly emitting a photon. 
This is a new class of dibaryons with decay widths < IkeV. The contribution of such dibaryons to strong interaction 
processes of hadrons is small. However, their contribution to electromagnetic processes on light nuclei may exceed the 
cross section for the process under study out of range of the dibaryon resonance by several orders of magnitude ||,0,D . 
I 1 The experimental discovery of such states would have important consequences for particle and nuclear physics. 

In the frame of the MIT bag model, Mulders et al. Q calculated the masses of different dibaryons, in particular, 
NN-decoupled dibaryons. They predicted dibaryons D(T = 0; J p = CT,1", 2~;M = 2.11 GeV) and D(l;l";2.2 
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GeV) corresponding to the forbidden NN states 13 Pj and 31 Pi. However, the dibaryon masses obtained exceed the 
■ pion production threshold. Therefore, these dibaryons can decay into the irNN channel. The possibility of the 
existence of dibaryons with masses M < 2mjv + was discussed by Kondratyuk et al. J^] in the model of stretched 
rotating bags, taking account of spin-orbital quark interactions. In the frame of the chiral soliton model, Kopeliovich 
[ fTpf predicted that the masses of D(T = 1, J p = 1 + ) and D(0, 2 + ) dibaryons exceeded the two-nucleon mass by 60 
and 90 MeV, respectively. These values are lower than the pion production threshold. 

Unfortunately, all results obtained for the dibaryon masses are model dependent. Therefore, only an experiment 
could answer the question about the existence of SNDs and their masses. 

In the work [ pd[ the existence of a dibaryon, called d' , with quantum numbers T = even and J p = 0~ which 
forbid its decay into two nucleons, and with the mass M = 2.06 GeV, and the decay width T^nn — 0.5 MeV, has 
been postulated to explain the observed resonance-like behavior in the energy dependence of the pionic double charge 
exchange on nuclei at an energy below the A-resonance. However, there is a more convential interpretation of these 
data. It was shown jl2| in the frame of the distorted-wave impulse approximation that such a peak arises nuturally 
because of the pion propogation in the sequential process, in which pion double charge exhange occurs through two 
successive nN charge exchange reactions on two neuterons. 

In the work |p| dibaryons with exotic quantum numbers were searched for in the process pp — > ppjj- The 
experiment was performed with a proton beam from the JINR phasotron at the energy of 198 MeV. Two photons 
were detected in the energy range between of 10 and 100 MeV. Some structure was observed in the photon energy 
spectrum, which was attributed to the exotic dibaryon with the mass M w 1920 MeV. The statistical significance of 
the effect is about 8c. However, a big width (FWHM=31.3±5.0 MeV) of the observed structure was obtained, and an 
additional careful experimental study of this reaction is needed in order to understand the nature of this structure. It 
is necessary to measure the photon energy spectrum with a higher energy resolution and a higher statistical precision. 
Moreover, measurements at two different incident proton energies are needed to check whether the structure observed 
is a SND. 
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On the other hand, an analysis fllif of the Uppsala proton-proton bremsstrahlung data looking for the presence of 
a dibaryon in the mass range from 1900 to 1960 MeV only gave the upper limits of 10 and 3 nb for the dibaryon 
production cross section at proton beam energies of 200 and 310 MeV, respectively. This result disagrees with the 
data 1 1 3 1 . 

Earlier we have carried out measurements of missing mass spectra in the reaction pd — > pX [El| using a double-arm 
spectrometer at the Moscow Meson Factory. The narrow peak in the missing mass spectrum at 1905 MeV with a 
width equal to the experimental resolution of 7 MeV has been observed in this experiment. 

As will be shown below, the nucleons and the deuteron from the decay of the dibaryons under consideration into 
jNN and jd have to be emitted in a narrow angular cone with respect to the direction of motion of the dibaryon. 
The size of this cone depends also on the quantum numbers of the dibaryon. So a detection of the scattered proton 
in coincidence with the proton (or the deuteron) from dibaryon decay at correlated angles gives a good possibility to 
separate the SNDs from the background and to determine their quantum numbers. 

In the present paper we give a more detailed study of the reactions pd — ► p + pX\ and pd — > p + dX^ (where X\ 
and Xi are undetected particles in this experiment) with the aim to search for SNDs at the Moscow Meson Factory 
using an improved facility. We will consider the following dibaryons: D(T — 0,J P — + ), D(0,0 - ), £>(1,1 + ), and 
D(l,l"). 

It is worth noting that the state (T = 1, J p = 1 ) corresponds to the states 31 Pi and 33 Pi in the NN channel. The 
former is forbidden and the latter is allowed for a two-nucleon state. In our work we will study the dibaryon D(l, 1 _ ), 
a decay of which into two nucleons is forbidden by the Pauli principle (i.e. 31 Pi state). 

The contents of the paper are the following. In Sec. |l| the decay widths of SNDs are given and the cross sections 
of these dibaryon production in the processes pd — ► pD are calculated. A description of the experimental setup is in 
Sec. [II. In Sec. IV the results of the measurements are presented. The results of a Monte Carlo calculation and an 
analysis of the experimental data obtained are presented in Sec. M. The main conclusions are given in Sec. VI. 



II. CROSS SECTIONS OF THE SUPERNARROW DIBARYON PRODUCTION IN THE REACTION 

pd — » pD 

In the process pd — > pD, SNDs can be produced only if the nucleons in the deuteron overlap sufficiently, such that 
a six-quark state with deuteron quantum numbers can be formed. In this case, an interaction of a photon or a meson 
with this state can change its quantum numbers so that a metastable state is formed. Therefore, the probability of 
the production of such dibaryons is proportional to the probability r\ of the six-quark state existing in the deuteron. 

The magnitude of r\ can be estimated from the deuteron form factor at large Q 2 (see, for example, [0). However, 
the values obtained depend strongly on the model of the form factor of the six-quark state over a broad region of Q 2 . 
Another way to estimate this parameter is to use the discrepancy between the theoretical and experimental values of 
the deuteron magnetic moment [[l7 18 1. This method is free from the restrictions quoted above and gives rj < 0.03 



We assume here that rj and the probability of the full overlap of nucleons in the virtual singlet state 31 So are 
equal to 0.01. 

Since the energy of nucleons, produced in the decay of the dibaryons under study with M < 2toat + m n , is small, 
it may be expected that the main contribution to a two-nucleon system should come from 13 Si (deuteron) and 31 So 
(virtual singlet) states. The results of calculations of the decay widths of the dibaryons into 7c? and jNN on the basis 
of such assumptions are listed in Table |. 

As shown in || the main contribution to the dibaryon photoproduction from the deuteron is given by the one-pion 
exchange diagram. Therefore, we will describe the production of SNDs in the process pd — > pD using the one meson 
exchange diagram also. 

The production of the isoscalar dibaryons D(0,0 + ) and D(0,0~) will be calculated with the help of the pole 
diagram with w-meson exchange between the incident proton and the deuteron. The vertices D(0, + ) — > Lud and 
D(0, 0~) — > uid can be written as 
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f (0,0+) 



= %Vve^x*G^°, (1) 
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r D(o,o-)^rf = j^VvG^i", (2) 

where G p „ = u M ?v — r^Vu, ^i^u = w^p v — p^w^. Here v and w are the four-vectors of the deuteron and u> meson 
polarization, r and p are their four-momenta. 
As a result of the calculations we have 
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R i+) = R^+tm^mif+t), 

t = 2(m 2 N - E 1 E 2 +pip 2 cosfl), 

Ei,pi(E2,P2) are the energy and momentum of the initial (final) proton, 9 is the angle of the scattered proton 
emission, nid is the deuteron mass, g^jvw i s the constant of the tuN interaction, f 2 2 /^ are the coupling constants in 
the vertices L>(0, 1+) + ui — >• L>(0, 0+) and D(0, 1+) + u -> £>(0, 0~ ), respectively, and D(0, 1+) is the six-quark state 
with quantum numbers of the deuteron. 

The isovector dibaryons D(l, 1 + ) and _D(1,1~) can be formed mainly as a result of pion exchange between the 
incident proton and the deuteron. The vertices D(l, 1 + ) — > 7r + <i and .D(l, 1 — ) — > 7r + ci are written as 

r D(i4+)^ 7r d = |^V^ e A^A ( rG' p,y $2' T : ( 5 ) 

F D(1,1-)^^ = JfVV G ^®3 U , ( 6 ) 

where = £1^1^ - qi^iv, = £2^2^ - <72^£2i/, and £ 1(2 ) and <j 1(2 ) are the four-vectors of polarization and 

four-momenta of the -D(l, 1 + ) (13(1, 1~)), respectively. 

The calculations yield the following expressions for the cross section of these dibaryon production 

= |U[(M 2 + m 2 d - *) 2 - 4m 2 M 2 ], (7) 
dcr pd ^ p D{i-) _ 1 gl „ U ,, T 2 
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and 3i2/(47'") are the coupling constants in the vertex for the transition of the six-quark state Z?(0, 1 + ) into the 
dibaryons under consideration via the pion exchange. 

For a numerical calculation of the dibaryon contribution let us assume that g 2 NN /Aii = 14.6, Q^nnI^ ~ 19.2, 
r\ = 0.01. The constants g 2 2 /47r and / 2 2 /47r are unknown. They are strong interaction constants. In order not to 
overestimate the contributions of the dibaryons we put these constants equal to 1. The results of the calculation of the 
SND production cross sections in the process pd — > pD at the kinetic energy of the incident proton T\ = E\ — m = 305 
MeV and the emission angle of the scattered proton 9 — 70° as a function of the dibaryon mass M are shown in 
Fig. |l|. Here the solid , dashed, dotted and dash-dotted lines correspond to D(0, + ), D(0, 0~), D(l, 1 + ), and D(l, 1~) 
dibaryon production, respectively. 

As indicated above, the SNDs decay mainly emitting a photon. Therefore, if we limit ourselves to the investigation 
of pd interactions with the photon in the final state, then the contribution of such dibaryons will essentially exceed the 
cross sections of the background processes. On the other hand, the special choice of the pd — » pX kinematics allows 
us to allocate the area where the contribution of the SNDs dominates, even without detection of the final photon. 
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III. EXPERIMENTAL SETUP 



The reaction p + d — > p + X was studied at the proton accelerator of the Moscow Meson Factory at 305 MeV. A 
proton beam with an average effective intensity ~ 0.1 nA bombarded alternatively CD2 and 12 C targets of 0.14 and 
0.18 g/cm 2 , respectively. The pd reaction contribution was determined by subtracting the 12 C spectrum from the 
CD2 spectrum. The exposition time of the experiment using CD2 as a target was 100 hours. This period consisted 
of two runs corresponding to two different kinematic conditions in order to avoid a systematic error. The charged 
particles produced in the reaction in question were detected at different correlated angles by the Two Arms Mass 
Spectrometer (TAMS). 

Our layout is shown schematically in Fig. ||. The left movable spectrometer arm, which is a single telescope 
AE — AE — E, was used to measure the energy and the time of flight of the scattered proton at a fixed emission 
angle 0^ = 0. A remote control drive allowed to place this arm at various angles between 65° and 85°. In the present 
experiment TAMS detected the scattered proton at the angle 0l — 72.5° (and 70° in another run) in coincidence with 
the second charged particle (either p or d) from the decay of the particle X during 55 (and 45) hours. 

The detection of the second charged particle in the right arm at angles close to the emission angle of particle X with 
mass M allows to suppress essentially the contribution of background processes and increase the relative contribution 
of a possible SND production. The right arm consisted of three telescopes which were located at Br = 33°, 35° and 
37°. These angles correspond to directions of the emission of the dibaryons with a certain mass. 

Each telescope included two thin plastic scintillation AE detectors (4 x 4 x 0.5 cm 3 ) and one E detector of thick 
plastic (5 x 5 x 20 cm 3 , used in the right arm) or BGO crystal (5x5x8 cm 3 , used in the left arm). Each detector 
was viewed by a PMT-143 phototube coupled to specially developed fast electronic modules. A trigger was generated 
by 4-fold coincidence of the two AE detector signals of the left arm combined with those of any telescope of the 
right arm. A time resolution better than 0.5 ns was achieved, which allowed to suppress essentially an accidental 
coincidence background. The scattered proton E signals, selected in the coincidence, formed the energy spectrum and 
accordingly the missing mass spectrum. Each useful event including two time of flights and two energies were stored 
event by event and then analyzed off line. 

The elastic pd scattering was measured at various angles of the spectrometer arms. The proton energy spectra 
obtained were used to calibrate the spectrometer in the proton energy and, accordingly, in the missing mass. The 
missing mass resolution of the spectrometer was AM « 3 MeV, the angular resolution was 1°. 



IV. RESULTS OF THE MEASUREMENTS 

The experimental missing mass spectra obtained with the deuteried polyethylene target are shown in Figs. ||(a,b,c). 
Each spectrum corresponds to a certain combination of outgoing angles of the scattered proton and the second charged 
particle. These combinations in Figs. ||(b,c) are consistent with the shift of the emission angle Or of the dibaryon with 
the given mass when the angle 0l changes from 70° to 72, 5°. As is evident from Figs. ||(a) and||(b), a resonance-like 
behavior of the spectra is observed for the CD2 target in two mass regions at 1905 ± 2 MeV and 1924 ± 2 MeV. On 
the other hand, the experiment with the carbon target resulted in a rather smooth spectra |l9| . This smoothness is 
caused by an essential increase of the contribution of background reactions in the interaction of the proton with the 
carbon. 

It is worth noting that the resonance structure at 1905 MeV is observed in both runs at 9l = 70° and 9l = 72.5° 
and at different setup modifications. This essentially decreases the possibility of a random origin of the observed 
structure. The results obtained at different 6l were summed up to improve statistics. 

The mass spectra obtained in Figs. p|a,b,c) were interpolated by second order polynomials (for the background) 
plus Gaussians (for the peaks in Figs.^(a) and ^b)). Let us regard Fig. |(b). In this case the interpolation gave 
X 2 /( n — 5) = 0.81. We determine the number of standard deviations (SD) as 



^N eff + N b ' 

where N e f f is the number of events above the background curve and Nb is the number of events below this curve. 
Taking 5 points for the peak at M = 1905 MeV, we have 65/V185 — 4.8 SD. That corresponds to the probability of 
statistic fluctuations P equal to 2.7 x 10~ 5 |20| ]. 

Using the same procedure for the peak at M = 1924 MeV in Fig. |^(a) results in 4.9 SD and P = 1.6 x 10~ 5 . 

If, for the number of SD, we use the expression [pT| 
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Er 5 (Nu-N bi )/a? 



where AT ti , A^j, and Ui correspond, respectively, to total, background, and error bar data, then we obtain 4.2 S.D. for 
the peak at M = 1905 MeV and 4.4 S.D. for the peak at M = 1924 MeV. 

The widths of both observed peaks correspond to the experimental resolution (3 MeV) . 

The peak at 1924 MeV was obtained only for one spectrum close to the upper limit of the missing mass. In the 
other cases this mass position was beyond the range of measurement. Therefore, in the present work we analyze in 
more detail the peak at 1905 MeV only. 

The experimental missing mass spectra in the range of 1872-1914 MeV, after subtracting the carbon contributions, 
are shown in Figs. f|(a,b,c). 

As seen from Figs. || and [|, the resonance behaviour of the cross section exhibits itself in a limited angular region. 

V. ANALYSIS OF THE EXPERIMENTAL DATA 

If the observed structure at M = 1905 MeV corresponds to a dibaryon decaying mainly into two nucleons, then 
the expected angular cone size of emitted nucleons would be about 50°. Moreover, the angular distributions of the 
emitted nucleons are expected to be very smooth in the angular region under consideration. Thus, even assuming that 
the dibaryon production cross section is equal to that of elastic scattering (40 /16/sr) its contribution to the missing 
mass spectra in Figs. ||(a,b,c) would be nearly the same and would not exceed 1-2 events. Hence, the peaks found 
are hardly interpreted as a manifestation of the formation and decay of such states. 

It was shown in || that the decay of the SND into ^NN had to be characterized in the rest frame by a narrow 
peak near the maximum photon energy in the probability distribution of the dibaryon decay over an emitted photon 
energy. It leads to an essential limitation of the outgoing nucleon angles. On the other hand, if such a dibaryon decays 
into jd, the emitted deuteron angles are limited by the following condition: sin#d < Mp* d / (rridPD), where pu is the 
momentum (in the lab. syst.) of the dibaryon and p* d is the momentum (in the c.m.s.) of the deuteron. 

Using the Monte Carlo simulation we estimated the contribution of the SNDs with different quantum numbers and 
M=1905 MeV to the mass spectra at various angles of the left and right arms of our setup. The production cross 
section and branching ratio of these states were taken according to the calculations presented in Sec. [D] and the proton 
beam current was assumed to be equal to 0.1 nA. The results obtained for different quantum numbers and decay 
modes of the dibaryons are listed in Table ^ 

This calculation showed that the angular cone of protons and deuterons emitted from a certain dibaryonic state can 
be rather narrow. The axis of this cone is lined up with the direction of the dibaryon emission. Therefore, placing the 
right spectrometer arm at the expected angle of the dibaryon emission we increase essentially the signal-to-background 
ratio. 

Figs. [5](a,b) exhibit the angular distributions of charged particles (either p or d) from the decay of the dibaryon 
D(l,l + ) (Fig. ||(a)) and D(0,0 + ) (Fig. |^(b)). The solid and dashed curves correspond to the emission angles of 
the scattered proton equal to 70° and 72.5°, respectively. The solid vertical lines show the location of the right 
arm detectors. These figures demonstrate the dependence of the angular distribution of the charged particles under 
consideration on the quantum numbers of the dibaryons. 

In Figs. |i](a,b,c), the experimental spectra are compared with the predicted yields normalized to the maximum 
of the measured signal in Fig. |I](b). The solid and dashed curves in these figures correspond to the states with the 
isospin T=l and T=0, respectively. The background is described by a second order polynomial. 

As is seen from these figures and Table O, the ratios of the calculated yields to the given spectra are expected to be 
0.3 : 1 : 0.7, if the state at 1905 MeV is interpreted as an isovector dibaryon (D(T = 1, J p = 1+) or D(l, 1")). This 
is in agreement with our experimental data within the errors. On the other hand, the signals from isoscalar dibaryons 
(D(0,0 + ) or D(0,0~)) could be observed in Figs. |J(b,c) with the same probability. So, the dibaryon with M = 1905 
MeV is most likely to have the isospin equal to 1. 

In order to estimate the dibaryon production cross section, we use normalization to the elastic pd scattering and 
then, comparing the theoretical predictions for the yields of one of the isovector dibaryons with the experimental data, 
find that the differential cross section of this dibaryon production is equal to 8 ± 4 pb/ sr. That allows us to evaluate 
the product of unknown constants in Eqs. (Q) and (g). If it is a dibaryon D{1, 1 + ), then r]g1/4ir = (0.44±0.22) x 10~ 2 . 
But if it is a dibaryon D(l, 1~), we have rjg^/^n = (0.2 ± 0.1) x 10~ 2 . The quantum numbers J p of the observed 
state could be found, for example, by a search for such dibaryons in the processes of charged pion photoproduction 
by polarized photons from the deuteron M . 
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As for the structure at 1924 MeV, the number of events in this peak exceeds the expected yields of the usual 
dibaryons (decaying into two nucleons) by almost two orders, even when produced with the rather big cross section 
of 40/i6/sr. So this state, probably, could be a supernarrow dibaryon, too. 

It should be noted that the reaction pd — > NX was investigated in other works, too (see, for example, p2|). 
However, in contrast to the present work, the authors of these works did not study either the correlation between 
the parameters of the scattered nucleon and the second detected particle or the emission of the photon from the 
dibaryon decay. Therefore, in these works the relative contribution of the dibaryons under consideration was small, 
which hampered their observation. 



VI. CONCLUSIONS 



The following conclusions can be made . 

1) As a result of the study of the reactions pd — > ppX\ and pd — ► pdX-i two narrow peaks at 1905 and 1924 MeV 
with widths less than 3 MeV have been observed in the missing mass spectra. 

2) The analysis of the angular distributions of the protons and the deuterons from decay of particle X, produced 
in the reaction pd — > pX, showed that the peak found at 1905 MeV can be explained as a manifestation of the SND 
with isospin equal to 1, the decay of which into two nucleons is forbidden by the Pauli exclusion principle. 

3) Probably, the observed state at 1924 MeV can be interpreted as the SND, too. 
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TABLE I. Decay widths of the dibaryons D(0, + ), D(0, 0"), D(l, 1+) and 15(1, 1") at various dibaryon masses M. 







D(0,0+) 






0(0,0") 




£>(1,1+) 


75(1,1") 


M 






r t 


r 7 d 


r^pri 


r t 


r t ~ r 7 jvjv 


r t ~ t 7 jvjv 


(GeV) 


(keV) 


(keV) 


(keV) 


(keV) 
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TABLE II. The expected contributions of the SNDs with M=1905 MeV and different quantum numbers to the spectra of 
mass at various angles of the left and right arms of the setup. 





6l 






70° 










72.5° 








e R 


31° 


33° 


35° 


37° 


39° 


31° 


33° 


35° 


37° 


39° 


T, J p 
























0,0+ 




5 


8 


8 


7 


5 


9 


10 


10 


6 


3 


0,0+ 


'yd 





38 


17 


32 





4 


23 


34 








0,0~ 


7pn 


3 


4 


5 


4 


2 


4 


6 


5 


2 


2 


0,0" 


'yd 





19 


10 


17 





2 


12 


17 








1,1+ 


jpn 


14 


55 


81 


50 


12 


47 


97 


79 


23 


5 


1,1" 


'ypn 


30 


120 


180 


110 


27 


103 


213 


174 


50 


12 



8 




1880 1900 1920 1940 1960 1980 2000 

M D , MeV 

FIG. 1. The cross section of the supernarrow dibaryon production in the process pd — > pD at Ti = 305 MeV and 9 = 70° 
as a function of the dibaryon mass M. The solid, dashed, doted and dash-doted lines correspond to the production of the 
dibaryons Z>(0,0+), D(0,0"), D(l, 1+) and D(l, 1"), respectively. 
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FIG. 2. The Two Arm Mass Spectrometer TAMS. So, Si, S2 and S3 are start detectors; Fo, Fi, F2 and F 3 are stop AE 
detectors; Do is a BGO detector; Di, D2 and D3 are full absorption E detectors. 
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FIG. 3. The missing mass spectra of the reaction on CD2: (a) from the run data at 6l = 70° and Or — 33°, b) from the run 
data at 70° and 35°, summarized with the run data at 72.5° and 33°, c) from the run data at 70° and 37°, summarized with 
the run data at 72.5° and 35°. 
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FIG. 4. The missing mass spectra of the reaction on deuteron. Areas marked by letters a), b), c) correspond to the same 
experimental conditions as in Fig. 3. The solid and dashed curves are normalized theoretical predictions of the yields of the 
supernarrow dibaryon with isotopic spin equal to 1 and 0, respectively. 
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FIG. 5. The angular distribution of protons emitted by the dibaryon D(l, 1 + ) (a), and the sum of the angular distributions 
of protons and deuterons emitted by the dibaryon D(0, + ) (b). The solid and dashed curves correspond to the scattered proton 
angles 70° and 72.5°, respectively. The solid vertical lines show the location of the right arm detectors. 
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